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METHOD OF DETERMINING THE IMPEDANCE OF AN ELECTROCHEMICAL SYSTEM 



To determine properties and qualities of an electro- 
chemical system (such as a cell,, a battery, an electrodeposition 
system, a system for analyzing a medium) and predict its future 
operation, one of the basic parameters is the system's 
5 impedance . 

Generally, an electrochemical system comprises, in an 
electrolytic medium, two main electrodes - a work electrode and 
a back- elect rode. A reference electrode is arranged in the 
vicinity of the work electrode and is used in relation therewith 

10 to perform various measurement or regulation operations. In the 
following, the mentioned voltages and currents correspond to 
measurements performed between a work electrode and the corre- 
sponding reference electrode. 

The impedance of an electrochemical system corresponds 

15 to the voltage/current ratio. It is known that to rapidly 
measure the value of this impedance according to frequency, it 
is desirable to calculate the ratio of the discrete Fourier 
transforms (DFT) of the voltage and of the current. For this 
purpose, an excitation signal of small amplitude is applied 

20 between a reference electrode and a work electrode of an 
electrochemical cell, and N successive values e(n) and i(n) of 
the voltage and of the current, with 1 < n < N, are measured at 



equal time intervals AT. The general expression of discrete 

Fourier transforms E(f m ) and I(f m ) of the voltage and of the 
current for a sequence of N points is: 

E(f m ) = AT Je(n)exp(-27cjf m nAT) (1) 



n=l 
N 



5 I(fJ = ATXi(n)exp(-2rcjf m nAT) (2) 

n=l 

i, 

where j designates the complex number having -1 as a square, N\ 
designates the number of measurement points, AT the sampling 
interval. The calculation of the discrete Fourier transform may 
be performed for N/2 values of frequency f m , with 0 < m < (N/2)- 
10 1, N being an even number. These N/2 discrete frequencies are 
regularly distributed between 0 and 1/2AT (0. . .1/mAT. . .1/2AT) . 

The complex impedance for a given frequency f m is then 

equal to: 

Z(f m ) = E(f m )/I(f m ). (3) 

15 The impedance measurement methods used in practice, 

mainly differ by the type of excitation signal of the system: 
sinusoidal, multi- sinusoidal, white noise, etc. A sinusoidal 
excitation is by far the most widely used, since it appears to 
be the most accurate. The method consists of imposing a 

20 sinusoidal regulation to the electrochemical cell with a signal 
of small amplitude and of recording the cell current and voltage 
response. The ratio of the Fourier transforms of the voltage and 
of the current at the sinus frequency will give the impedance 
value at this frequency. The frequency spectrum is swept by 

25 modifying the excitation frequency- 

In the theory of electric systems, the expression of 
the impedance (equation (3)) is correct, provided that the 
analyzed system is linear and stationary and that disturbances 
are not introduced by phenomena external to the system. In the 

3 0 case of electrochemical systems, the fulfilling of these condi- 
tions imposes specific precautions. The transfer function of an 
electrochemical system is generally non- linear but it can be 
considered as linear over a small portion. This is why a low- 
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amplitude excitation signal is used. Accordingly, the signal-to- 
noise ratio decreases and the measurement time, that is, the 
number of measurement points, must be increased to suppress the 
noise by an integration of the response. However, the 
5 lengthening of the measurement time sets problems regarding to 
the stationarity of the electrochemical system. In many real 
cases, the stationarity condition cannot be fulfilled. The 
causes are many: voltage relaxation, current relaxation, 
concentration relaxation. In such conditions, if the measurement 

10 time is sufficiently long to have a good signal-to-noise ratio, 
the system is not stationary during the measurement time and the 
calculated impedance has no great meaning any more, especially 
at low frequencies (for example, smaller than one hertz) where 
non- stationarity problems are particularly serious. 

15 Various methods have been provided to solve this 

problem of the non- stationarity of electrochemical systems, but 
none has yielded a satisfactory solution. 

Thus, an object of the present invention is to provide 
a novel method for calculating the impedance of an electro- 

20 chemical system enabling ignoring errors linked to the system 
non- stat ionarity . 

To achieve this object, the present invention provides 
a method for determining the complex impedance Z(f m ) of a non- 
stationary electrochemical system, comprising the steps of: 

25 setting the system to a selected voltage state and 

applying a sinusoidal signal of frequency f m thereto, 

measuring, immediately after, successive values of the 
voltage and of the current at regular time intervals AT, 

calculating the discrete Fourier transforms of the 

3 0 voltage (E(f)) and of the current (1(f)), the voltage transform 
being calculated for the sole frequency f m of the sinusoidal 
signal and the current transform being calculated for frequency 
f m and for two adjacent frequencies f m -i and f m+ i on either side 
of frequency f m , and 
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calculating the impedance according to the following 

formula: 

Z(f m )=E(f m )/l*(f m ) 
where I* designates a corrected current such that: 
5 Re[I*(f m )] = Re[I(f m )] - {Re [I (f m+1 ) ] +Re [I (f^) ] }/2 

Im[I*(f m )] = Im[I(f m )] - {imfKf^xJl+ImtKf^!)] }/2. 
A specific embodiment of the present invention will be 
non- limit ingly discussed in relation with the appended drawing 
which shows the amplitude spectrum according to the frequency of 
10 an electrochemical system submitted to a sinusoidal excitation 
of small amplitude and to a voltage step. 

It should first be noted that, since the applied volt- 
age is imposed by the analysis tool, the system non-stationarity 
can only appear as current variations. 
15 The present invention is based on the analysis of the 

behavior of an electrochemical system submitted to a voltage 
step. The DFT of the current response mainly translates relaxa- 
tion phenomena and thus the effect of the system non- 
stationarity. As shown by the curve in dotted lines on the 
20 single drawing, the spectrum of amplitude I of the unit of the 
current DFT usually exhibits a strong low- frequency response to 
voltage steps. 

Conversely, the response to a sinusoidal excitation of 
frequency f m of a stationary system translates as a single line 
25 at frequency f m . In practice, this response corresponds on the 
one hand to the system response to the excitation at frequency 
f m , and on the other hand to the contribution of relaxation 
effects. 

The advantage of the vise of the Fourier transform is 
3 0 that the single -frequency spectrum linked to the sinusoidal 
excitation is superposed to the spectrum linked to the voltage 
step. 

According to the present invention, once the above 
considerations have been taken into account, it is provided to 
3 5 subtract from the intensity response at frequency f m the contri- 
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bution due to relaxations, evaluated based on the analysis of 
the DFT at frequencies adjacent to f m . Indeed, the Fourier 
transform for frequencies f m -i and f m+1 of the system will only 
correspond to non-stationarity phenomena and it will be 
5 considered that the amplitude value linked to non-stationarities 
for a frequency f m is the average of the values for the two 
adjacent frequencies. 

Thus, the present invention provides: 

- setting an electrochemical system to a selected 
10 voltage state and applying a single -frequency sinusoidal excita- 
tion thereto, 

- measuring, without waiting for the system stabili- 
zation, the amplitude of the current and of the voltage at 
regular intervals immediately after application of the voltage 

15 while the sinusoidal excitation is applied, and 

- calculating, on the one hand, the discrete Fourier 
transform of the voltage for frequency f m of the excitation and 
on the other hand the DFT of the current for value f m and for 
two frequencies f m -i and f m +i adjacent to frequency f m , and 

20 - calculating the value of the complex impedance at 

frequency f m based on the value of the TFD of the voltage for 
frequency f m and on a corrected value of the TFD of the current, 
taking into account the TFDs calculated for frequencies f m -i# 
f m , and f m+1 . 

25 Corrected amplitude I* of the discrete Fourier trans- 

form of the current will be calculated in real value and in 
imaginary value by the two following equations: 

Re[I*(f m )] =Re[I(f m )] - {Re[I(f m+1 )+RetKf m -i)]}/2 
Im[I*(f m )] = Im[I(f m )] - {lm[I(f m+1 )^Jm[I(f m . 1 )]}/2. 
30 The value of the impedance corrected with the non- 

stationarity effects then is: 

Z(f m )=E(f m )/l*(f m ) 
with I*(f m ) = Re[I*(f m )]+jI m [I*(f m )] 



As soon as the reading of the points is performed for 
frequency f m , a sinusoidal signal can be applied to a new 
frequency and a new reading can be performed, and so on. 

An advantage of the present method is that it makes a 
correct impedance analysis on non- stationary electrochemical 
systems possible, especially at very low frequencies. At the 
same time, a considerable time gain is obtained for systems 
which stabilize slowly since with the correction according to 
the present invention, it is no longer necessary to wait for the 
stabilization after powering-on of the system to start an imped- 
ance analysis. 



